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INTRODUCTION 

This is an introductory course on the open problems of fully developed turbulence 
which present a long standing challenge for theoretical and mathematical physics. The 
plan of the course is as follows: 

Lecture 1. Hydrodynamical equations. Existence of solutions. Statistical description. 
Kolmogorov scaling theory. 

Lecture 2. Functional approach to turbulence, similarities and differences with field 
theory. 

Lecture 3. Passive scalar and breakdown of the Kolmogorov theory. 
Lecture 4. Inverse renormalization group. 

LECTURE 1 

The hydrodynamical flows in gasses and liquids are believed to be described in a 
variety of realistic situations by incompressible Navier-Stokes (NS) equations 

dtv + {vV)v -uAv = -if -Vp), V-v = (1) 



where the vector v{t, x) describes the velocity of the fluid at time t and space point x, 
positive constant z/ is the viscosity, p is the constant density^, f(t,x) is the (intensive) 
force and p(t, x) is the pressure. The above equations date back to the works of Navier 
(1823) and Stokes (1843) and modify the even older Euler equation (1755) by the 
addition of the dissipative term uAv. In most physical applications the dimensionality 
of the space is 3 or 2, but the equations make sense in general dimension d. The Euler 
equation without the uAv term is really the F = ma (or rather a = -p;^F) equality 
for the volume element of the fluid. The uAv term represents the friction forces. By 
taking the scalar product of the NS equation with v and assuming that the flow 
velocity vanishes sufficiently fast at large distances (or at the boundary), one deduces 
the energy balance: 

dt\jv^ = -vj{Vvf + Jf-v. (2) 

-^lectures given at the 1996 Cargese Summer Institute, July 22 - August 3 
^we shall set it to 1 below 



The equation says that the rate of change of fluid energy is equal to the energy injection 
J f - v (the work of external forces per unit time) minus the energy dissipation per unit 
time z/ /(Vf)^ due to the fluid friction. 

The Euler equation has a nice infinite-dimensional geometric interpretation: it de- 
scribes the geodesic flow on the group of volume preserving diffeomorphisms^ . We 
shall briefly sketch this argument. Let, more generally, the space be an oriented 
(i-dimensional Riemannian manifold {M,g). Denote by r the Riemannian volume: 
r(x) = \J g{x) dx^ A ■ ■ ■ A dx'^ . Let Dif fr be the group of diffeomorphisms of M 
preserving r: (f)*T = t. The space Vectr of divergenceless vector fields v (i.e. the ones 
preserving r, C^r = 0) may be considered as the Lie algebra of Diffr^. Upon the 
identification of infinitesimal variations of diffeomorphisms (p with the vector fields v 
by 5</)(x) = v{(t>{x)), the scalar product of vector fields 

IK'T = J 9iv,v)r = J gij{x)v\x)v^{x) t{x) 

M M 

induces on Dif fj. a right-invariant Riemannian metric. This metric in not left-invariant. 
Indeed, 

dL^dR^-i V = (j)^v 

where dL^ and dR^ stand for the tangent maps to the left and right translations by 
on Dif fr and 

(f)*v{x) = d(j)v{(j)~^{x)) 
is the pushforward of v by 0. But, in general, 

ii0*^ir ^ ii^ir 

since the tangent map dcj) does not preserve the Riemannian length of vectors. 

The geodesic flows with respect to left-right-invariant metrics on a group are, mod- 
ulo time-independent left and right translations, one-parameter subgroups. This is not 
the case if the metric is only right-invariant. The geodesic motions t ^ (f)t on Dif f^ 
corresponding to the right-invariant Riemannian metric defined above extremize the 
action 

^(0) = I / wvm'dt 

where 

dt(t){t,x) = v(t,(f){t,x)) = v(t,y). 

ExphcitlyQ, 

SS{cf>) = \5 j g,,{y)v\t,y)v^{t,y) r(y) dt 
■^recall that the Euler top describes the geodesic flow on the group SO{3) 

"^for non-compact M we shall assume that (p do not move points outside a compact subset and that 
V have compact support 

Afield theoriests will notice a similarity of the following calculation to those in nonlinear sigma 
models 



J gij{(t){t,x)) dt(j)'{t,x) dt(j)^{t,x) t{x) dt 
\ j dkgiji(pit,x)) S(p^{t,x) dt(p\t,x) dt(p^{t,x) t{x) dt 

+ J gij{(j){t,x)) dt(t)\t,x) dt5(jy'{t,x) r(x) dt 
JS(j)\t,x) \djgik{y)v\t,y)v^{t,y) - dkgij{y)v\t,y)v^{t,y) 



-g,,{y)[dtv\t,y) + v\t,y)dkv\t,y) 



T{y) dt 



9ij 



dtv' + v'^dkv' + W^'idkgin + digun - dngik)v''v^ 



T dt 



- j g{u, dt + VvV ) T dt 



where u^{t,y) = 6(j>'{t,x) is an arbitrary divergenceless vector field. 

{V^vy = w'dkv'' + {^^}w''v' 

denotes the covariant derivative of the vector field v in the direction of w with the 
Levi-Civita symbols 

{fcj = \g'''{dkgin + digkn- dngik) ■ 

Since a vector field on M orthogonal to all vector fields without divergence is a gradient, 
we obtain the (generalized) Euler equation 

dtv + VvV = — Vp 

for the divergenceless vector fields v. The NS term uAv still makes sense in the geo- 
metric setup and may be added to the equation. 

The Euler and NS equations are examples of nonlinear partial differential evolution 
equations. After a century and a half study, they still pose major open problems as far 
as the control of their solutions is concerned. The most interesting questions concern 
the short-distance (ultra-violet) behavior. Suppose that we start from smooth initial 
data and the force / is smooth. For simplicity, let us assume compact support of both 
(we may also consider the compact space or the case with boundary conditions). It 
is known that the smooth solutions of the so posed problems are unique and exist for 
short time. Do they exist for all times? The answer is positive in 2 dimensions for 
both Euler and NS equations [43, pH] but in 3 dimensions the answer is not known. 



It is expected to be positive in the NS case. The opinions about the Euler case (no 
blowup versus finite-time blowup for special smooth initial conditions) are divided and 
fluctuate in time. 

In an important 1933 paper on the NS equation Leray has introduced the notion 
of weak solutions of the equation. A vector field v{t, x) locally in is a weak solution 
if 



and 







for any smooth divergenceless vector field u and any smooth function y9, both with 
compact supports. Leray showed by compactness arguments existence of weak global 
solutions of the 3-dimensional NS equations with additional properties (e.g. with space 



derivatives locally square integrable). The weak solutions are not unique (there are 
weak solutions of the 2- dimensional Euler equation with compact support W^)- 

The NS equation is invariant under rescalings. Let 

v{t,x) = TS"^ v{Tt, Sx) , 
f {t,x) = t'^S'^ f{Tt,Sx), 
p{t,x) = t'^S"'^ v{Tt, Sx) , 
U = TS~'^ V . 

If V and p solve the NS equation with viscosity v and force / then v and p give a solution 
for viscosity v and force /. It is convenient to introduce the dimensionless version of 
the (inverse) viscosity, the Reynolds number 

V 

where b^y is a characteristic velocity difference over scale L of the order of the size of 
the system. For example, for the flow in a pipe of radius L, we may take b^v as velocity 
in the middle of the pipe (minus the vanishing velocity on the wall of the pipe). A basic 
phenomenological observation in hydrodynamics is that for R<\ the flows are regular 
(laminar) whereas i? ^ 1 correspond to very irregular (turbulent) flows with a rich set 
of scenarios occurring for intermediate R. Note the scale-dependent character of the 
Reynold number. Following define the running Reynolds number 

where 0{r) = { (s, y) | |s — t| < ^ , |?/ — x| < r } . Note that we may rewrite 

brV r 
z/ 

where SrV, the mean velocity difference on scale r, is calculated by multiplying the 
mean square gradient of v over 0{r) by r. The best regularity result about the weak 
solutions of the NS equation is due to Caffarelli-Kohn-Nirenberg |^ and says that there 
exist e > such that if Rr < e then the solution is smooth on the 0{er) neighborhood 
of (t, x) . This implies that for a weak solution, the Haussdorf dimension of the set 
of singularities is < 1. Note the spirit of the result in line with the phenomenological 
characterization of laminar flows. 

For high Reynolds numbers it is reasonable to attempt a statistical description 
of complicated turbulent flows. In the theoretical approach, the statistics may be 
generated by considering random initial data or/and random forcing. v{t,x) becomes 
then a random field. We shall be interested in describing a stationary statistical state 
of the latter. In such a state the mean overall energy of the fluid is constant in time so 
that the energy balance equation (||) implies that 

jiH^vr) = jif-v), 

where ( — ) denotes the ensemble average, or, in a homogeneous state. 



(3) 



where e denotes the mean dissipation rate and ^ the mean injection rate of energy, 
both with dimension '^^^3 ■ In the situation where the energy injection is a large 
distance process (e.g. in the atmospheric turbulence) one expects that for high R 
a scale separation occurs with the energy dissipation taking place on much smaller 
distances. Pictorially, energy is transmitted to the fluid by induction of large eddies on 
scale L which subsequently break to smaller scale eddies and so on passing energy to 
shorter and shorter scales without substantial loss until the viscous scale 77 is reached 
where the friction dissipates energy. Such a (Richardson [^) energy cascade is then 
characterized by the integral scale L, the viscous scale ri and the energy dissipation 
rate e. The scale ratio Ljri should grow with the Reynolds number. The interval of 
distance scales r satisfying L ^ r ^ is called the inertial range. The cascade picture 
may be formulated in more quantitative terms by introducing the quantities [1^ 



1<K = V J (^J{Vv{0)-Vv{x))e-'''-''dx^dk 



\k\<K 

(2^)- 



(dk = j^^) and 



= / (/ (/(O) ■ v{x)) e-''-^ dxyk 



\k\<K 

interpreted as the mean dissipation and mean injection rate in wavenumbers k with 
\k\ < K. The injection of energy limited to distances >L means that, as a function of 
ip<_K is close to e everywhere except for K<j^ where it falls to zero with K ^ {]. 
Similarly, the cascade picture should imply that e<x is negligable for <^ ^ and then 
grows to e. The difference 

has the interpretation of the energy flux out of the wavenumbers k with \k\ < K. It 
should be approximately constant and equal to e in the inertial range <^ K <^ ^ . 

In 1941, Kolmogorov |23| went one step further by postulating that the equal-time 



correlators of velocity differences over distances in the inertial range should be universal 
functions of the latter and of the dissipation rate e. In particular for the structure 
functions 

Sn{x) = {{v{x)-vm-xr) 

with 2=1^, assuming also isotropy of the turbulent state, one obtains 

with r = |x| and universal constants C„. Indeed, the right hand side is the only 
function of e and r with the dimension Cume )^- Kolmogorov theory implies that the 
typical velocity difference over distance r behaves as gi/^^i/s^ Pqj^ ^^j-^g scale-dependent 
Reynolds number one obtains then Rj. ~ e^^lrl!l_ particular, R = Rl ~ ^^^^l*^^ 
0{1) = Rrj ^ ^^^^J*^'' hence £ ~ i?"^/^ and it decreases with R. 

For n = 3 the relation (^) is essentially a rigorous result with C3 = — | in 3 
dimensions. Let us sketch the latter argument ||2^, |TB|. Suppose for simplicity that the 
force / is a random Gaussian field with mean zero and the covariance 

{f{t,x)f{s,y)) = 6{t-s)C^{x-y) (5) 



with diC"^ = . Inferring from the NS equation that 

t+At 

v(t + At) = v{t) - {{v- Vv)v - uAv + Vp)\,At + / /(s) ds + 0{Atf 

•I t 

the stationarity of {v(t,x) ■ v(t,y)) imphes that 

2z/ ( Vt'(x) ■ Vv{y)) = — z/( At>(x) ■ v{y)) — (x) ■ Av{y)) 

= — ((f (x) ■ V)f (x) ■ v{y)) — (f (x) ■ (f (y) ■ V)f (y)) + tr C(x — y) 

(the pressure does not contribute due to the transversahty of v). The first two terms 
on the right hand side may be rewritten as 

'^V,-{iv{x)-viy)fiv{x)-v{y))) 

so that we obtain the relation 

z/(Vw(x) ■ Vt;(0)) - jV^- {{v{x) - v{0)f{v{x) - v{0))) = itrC(x). (6) 

Taking first the hmit x ^ for positive u and assuming that the presence of the latter 
smoothes the behavior of ((f (x) — v{0)y{v{x) — f (0))) so that the second term on the 
left hand side vanishes when x — > 0, we obtain 

e = itrC(O) 

which is nothing else then the energy balance equation (^). On the other hand, taking 
first the u ^ limit of eq. (P) for x 7^ 0, we obtain 

- \V.-{{v{x)-v{0)nv{x)-vm)U = |trC(x). (7) 

The assumption that the force acts only on distances >L means that for r <^ L, 
C(x) = C(0) so that eq. (j^) means that in the inertial range 



\v , ■ {{v{x) - vm\v{x) - v{m 



e. 



Together with isotropy, this implies the relation 16 



{{v%x) - v\0)){v^{x) - v^{0)){v\x) -v''{0))) = -^(^{S'^x'' + + (5^V) 



from which the n = 3 case of eq. (1) follows with C3 



12 



d{d+2) • 



The structure functions are measured, more or less directly, in atmospheric or ocean 
flows, in water jets, in aerodynamic tunnels or in subtle experiments with helium gas 
inbetween rotating cylinders. They are also accessible in numerical simulations. One 
extracts then the scaling exponents assuming the behavior 

Sn{x) ~ r^" ■ 

(3 agrees well with the theoretical prediction (^3 = 1. Here are some other exponents 
obtained from wind tunnel data 

C2 = .70(.67), (4 = 1.28(1.33), (5 = 1-53(1.67), (6 = 1-77(2), (7 = 2.01(2.33) 



with the Kolmogorov values in the parenthesis for comparison. The discrepancy is 
quite pronounced (its direction is determined by the Holder inequality). One of the 
main open problems in the theory of fully developed turbulence is to explain, starting 
from the first principles (i.e. from the NS equation), the breakdown of the Kolmogorov 
theory leading to the anomalous structure-function exponents which indicate that the 
distribution of v{t) in the inertial range is far from Gaussian. The intuitive idea that 
only a part of fluid modes (temporal or/and spacial) participates in the turbulent 
cascade ("intermittency") has led to multiple models of the cascade of the (multi)fractal 
type . Such models are not based on the NS equation and allow to obtain essentially 
arbitrary spectra of (multifractal [|12|) exponents hence they do not really explain the 
breakdown of the normal scaling in realistic flows. 

For the "energy spectrum" 

\k\<K 

the Kolmogorov theory predicts 

e(ir)~ £2/3^-5/3 

for <^ K ^. The experimental data seem to confirm this behavior (with the 
possibility of a slight discrepancy consistent with the above value of ^2)- Deep in the 
dissipative regime {K ^ i), e{K) falls off much faster than in the inertial regime. 

LECTURE 2 

The NS equation with random force is an example of a stochastic evolution equation 
of the form 

9*$ = G(<l>) + / (8) 

with $ a function of time, G($) a functional local in time and / a stationary Gaussian 
process with mean zero and covariance 

(fit) f is)) = C{t-s). 

Upon solving the equation for given / with initial data $(to) = '^'o; one obtains a 
stochastic process t i— > $(t). The limit to ~^ —00, if exists, should then allow to pass 
to the stationary regime without dependence on the initial data. There is a simple way 
to write the expectation values of a functional of $ in terms of a formal functional 
integral: 

(^($)) = /^($) <5(a,$ - G($) - /) det(9,-§)D$rf/i,(/) (9) 

where dfi^ denotes the Gaussian measure with covariance C. Indeed, the integration 
over $ imposes the dependence (H) between $ and / and the / integral calculates the 
Gaussian expectation of JF($) viewed as a functional of /. The determinant in (||) is 
calculable (formally). 



det(9, - g) = det(9,) det(l-ar^f) 



1 xr- - y tr(9rl|g)" 

With the choice d^^{ti,t2) = 9{ti — t2), appropriate for solutions of the initial value 
problems, we have, setting^ 9{0) = i. 



tr 



5<S> 



for n > 1 

2-1 ^'^ 5'S> 



^ JtT^dt for n=l, 



where on the right hand side, |f is calculated at fixed time. Inserting the value of 
the determinant into eq. and rewriting the 5-function as an oscillatory integral, we 
obtain 

(^($)) = y'^($)e-(«.9*'f-G('f)-/)-|/tr|irf* DRD<l>dfiM) /norm., 
where (■, ■) stands for the scalar product, or, after integration over /, 

= J jr($) Q-siR,^) DJiD^ j ^orm. (10) 

where 

5(i?,$) = %{R, dt^-G{^)) + i 1 trgrft + |(i?,Ci?). (11) 

The above functional integral representation is known in the physical literature as the 
Martin-Siggia-Rose (MSR) formalism [Q. The field R is called the response field (its 
correlations measure the reaction of the system to small deterministic variations of the 
force). It may be integrated out in ([T0| ) leading to a Fokker-Planck type of functional 
integral but it will be more convenient to keep it in the functional representation. 

As a simple but instructive example, let us consider a linear stochastic equation 
describing forced diffusion: 

dtT = DAT + f (12) 

with the diffusion constant D and the forcing covariance 

{fit,x)f{s,y)) = 5{t-s)Cix-y) (13) 

where C is a smooth, fast decaying, positive-definite function. Given the initial data 
T(to), we may solve eq. (|1^): 

t 

T{t,x) = |e(*-*«)^^(x,y)T(to,t/)^^t/ + Jdsje^'-'^''^{x,y)f{s,y)dy 

to 

obtaining a Gaussian stochastic process with mean 

e(*-*o)^^T(to) 

and covariance 

((T(ti,Xi)-(T(ti,X i))) (T(t2, X2) - {T{h, X2))) ) 

^in simplest situations, this choice cancels in final expressions, in more complicated ones it is related 
to a choice of the (ordering) convention for the stochastic integrals 



min{ti,t2) 



to 

The limit to ~^ ~c>o results in a stationary Gaussian process with mean zero and 
covariance 

(T(ti,Xi)r(t2,X2)) = / ^\t.-t.\Dk--^H^.-^.) (14) 

Note that for ti = t2 the right hand side is a propagator of the scalar massless free field 
with UV cutoff C(A;). 

In the MSR formalism, the stochastic equation (|12D leads to the quadratic action 

S{R,T) = ijR{t,x){dtT - DAT){t,x)dtdx + \j R{t,x)C{x - y)R{t,y) dtdxdy 

which gives rise to the propagators 

{RR) = 0, 

{T{ti,xi)R{t2,X2)) = tidt + DA)-\h,xi;t2,X2) 

= ze(ti-t2)e(*^-*^)^^(xi,a;2) (15) 

and (TT) as in eq. ([I^). 

There are two interesting limiting regimes of forced diffusion: the one of long times 
and long distances and that of short times and short distances. To study the first one, 
let us introduce the rescaled fields 

Tx{t,x) = \^''-^^/^T{\t\\x), Rx{t,x) = \^''-^^^/^R{\H,\x) (16) 

for which 

{Uh,x,)Ut2,X2)) = X'-' j Ig^dk 

= f Q\ti~t2\Dk^-ik-(xi-X2) C(k/X) ^ £(0)_ f ^-\ti^t2\Dk^ + ikixi-X2) J_^j^ /^Y) 

J 2Dk X^oo 2-D y A; ' 

{TxRx) = {TR). 

The limiting (TxTx) covariance coincides with the one of the Langevin dynamics (of 
type A [^]) for the scalar massless field. Recall that the Langevin dynamics is given 
by the stochastic evolution equation 

with the delta-correlated noise 

{f{t,x)f{s,y)) = T5it~s)5{x-y). 

For V($) = J V{^(t)) dt, it describes convergence to the equilibrium (equal time) 
probability measure ~ e~^^'^^ . In our case 



v{T) = £^JiyT)\ r = c(o). 



The limiting behavior ([T7|) expresses the simple fact that at large distances the noise 
covariance C{x — y) looks like (/ C) 5{x — y). 

The study of the long-time, long-distance asymptotics is a typical field-theoretic 
problem where we are interested in the behavior of the system at distances much longer 
than the UV cutoff scale. Studying the opposite regime of distances much shorter than 
the cutoff scale may seem without interest. This is not so. In problems related to 
turbulence, function C{k) describes the spectrum of forcing rather than the momentum 
cutoff and we are interested in the regime of short times and short distances, more 
exactly in distance scales much shorter than the scale of energy injection interpreted in 
the field-theoretic context as the UV cutoff scale. Hence field theory and turbulence 
are concerned with opposite limiting regimes. In order to examine the short- 
time, short- distance asymptotics of forced diffusion, we rescale the fields differently 
introducing 

T^{t,x) = AT(^/A^a;/A), R^{t,x) = X''^-^ R{t / , x / X) . (18) 
and obtaining 

{T\tux,)T\t2,X2)) = X^J el*i-*^l^-'^^'-*'=(-i-2)^-^ i§,dk 

= ^'jl§h^^ - \m {\ti - h\ + i-,\x, - x,\') + 0{X') , (19) 

{T^R^) = {T R) 

for large A. Hence the (T^T^) covariance reaches a scaling form modulo a divergent 
constant. In other words, the differences T{t,x) — T{t,y) exhibit a scaling behavior at 
short times and short distances. The presence of the blowing up constant mode which 
may be eliminated by considering field differences is typical for the turbulence related 
problems. 

Returning to the NS equation (|^) with random Gaussian force (^ and applying the 
MSR formalism, we obtain the MSR action 

S{R,v) = i J R{t,x) ■ {dtv + v- Wv - i^Av){t,x) dtdx 
+ R{t,x) ■ C{x — y)R{t,y) dtdxdy 

with vector fields R and v satisfying V • i? = and V ■ f = 0. One could then set up a 
perturbative scheme by separating 

S{R,v) = So{R,v) + Si{R,v) 

with 

So{R, v) = ij R(t, x) ■ {dtv — z/At>)(t, x)dtdx + R{t, x) ■ C{x — y)R{t, y) dtdxdy^ 
Si{R,v) =ijRP^Xv-^v)dtdx, 

(P^j. stands for the projection on the transverse vector fields) and developing e""^^ in 
the power series, with 5*0 giving rise to the propagators 

{Rv), , {vv), 



and 5*1 to the vertex 



and with the simplest Feynman diagrams 



The perturbative expansion is plagued by the UV and IR divergences when v ^ and 
C — > const. Attempts were made to improve the situation by applying various resum- 
mation schemes, inspired by the field theoretic techniques (Schwinger-Dyson equations, 
renormalization group, etc.) but it seems fair to say that they were not very successful. 

Despite a formal similarity with the field theory formulation using functional inte- 
grals, the problem of the NS turbulence differs radically from the (dynamical formula- 
tion) of field theory in at least two crucial aspects: 

I. the nonlinear term v ■ Vf in the NS equation is not of the gradient type 
unlike the nonlinearities in the Langevin dynamics for field theory models, 

II. the noise (force) correlation in the NS turbulence should be close to a delta- 
function in the Fourier space rather than in the position space as for the Langevin 
dynamics. 

Consequently, the stationary statistical state of the NS fiuid is not of an equilibrium 
type unlike the Gibbs states corresponding to the stationary states of the Langevin 
dynamics for (euclidean) field theories. We have seen on the example of linear forced 
diffusion that the 2"'^ difference is already enough to make the problem very different 
from the field theoretic one. Let us profit from the occasion to mention that the study 
of the nonlinear version of forced diffusion: 

dtT = DAT - ^ + / (20) 

2C(0) ST J ^ ^ 

with random Gaussian force as in (|1^) is a problem of its own interest. The long- 
time long-distance asymptotics of the correlations is then a dynamical version of an 
interacting field theory problem (e.g. for V(T) ~ /T^) and may be studied by the 
renormalization group techniques. On the other hand, the short-time, short- distance 
asymptotics e.g. for V(T) ~ /(VT)'* is a completely open problemQ. The inverse 
renormalization group which will be the subject of the last lecture may provide a tool 
for the latter type of questions. 



LECTURE 3 

In view of the reputed difficulty of the NS problem, it would be instructive to 
consider simpler models of the stochastic evolution equation (§) randomly forced at 
long distances and with G($) not of the gradient type. Such a model is provided by 
an equation, with a relatively long history [^, [|, |T^, describing the passive advection 
of a scalar quantity T(t, x) (temperature) by a random velocity field 

dtT = -vVT + kAT + / . (21) 

^ we consider functionals of VT to avoid coupling to the constant mode unstable at short distances 
already in the linear case 



The positive coefficient k is called the molecular diffusivity. For a divergence-free w, 
operator — f ■ V is skew-symmetric (conserving the energy ^ / T^) and the —v-VT term 
is not of the gradient type, unlike the kAT one corresponding to a negative (energy 
dissipating) symmetric operator. Ideally, the velocity field v{t,x) should describe a 
turbulent NS flow but, to simplify radically the problem, one replaces it by a Gaussian 
random field. In Kraichnan noticed that in the case when v is decorrelated in time 



the problem becomes exactly soluble in the sense that one may write a closed system 
of differential equations for the correlation functions of T. In recent years Kraichnan's 
model of passive scalar has attracted a lot of attention, see e.g. [^, ^, |T^, as 
a prototype of a turbulent system in which one may study analytically the breakdown 
of a Kolmogorov-type scaling. 



We shall describe some results of the approach developed in |T8|, |19|, |^. The Gaussian 
velocity field v will be taken of mean zero and covariance 

{v\t,x)v^{s,y)) = 6{s-t)D'^{x-y) (22) 

with diD^^(x) = and 

D'^{x) = 2Dq5'^ - 2d'^{x) (23) 

where d^^{x) ~ for small r=|a;|, 0<^<2. Note the implied growth of the 2"^^ 
order velocity structure function with the distance, as in the turbulent flows. ^ = | 
corresponds to the Kolmogorov scaling dimension^ of v equal to |. More concretely, 
we may pose 



with a small IR regulator m, which leads to Dq ~ m^^ and 

d^^ix) ^ d^ix) - £j ((rf - 1 + OS'^ - r« = di^ix) (25) 

for small r and some constant D. Unlike Dq which diverges when m — 0, d]l^{x) 
possesses the m — limit equal to c?o (x) and scaling exactly with power ^. 

The source / will be assumed independent of v, Gaussian, with mean zero and with 
covariance 

{fit,x)f{s,y)) = 5{t-s)CC-^) 

for some positive definite, rotationally invariant test function C. 

The passive scalar model as set up above exhibits scale separation with the energy 
cascade in the inertial range of distances L > r > r/ = 0{{%f/^)L. The even structure 
functions of the scalar^ S'2„,(r) = {{T{t, x) ~T{t, 0))^") show exponential scaling in the 
inertial range. At least for small ^, 

^2n(r) ~ r^'" (26) 

with 

C2n = n{2-0-'-^^ + O{e)- (27) 



^as we shall see in the next lecture, the dimension of (5(i — s) is ^ — 2 
^the odd ones vanish 



In particular, the scaling of the 4**^ and higher structure functions is anomalous: the 
exponents (2n deviate from the ones of the Kolmogorov(-Corrsin) scaling theory giving 
C2n = n(2- Similar results were obtained in p, |^ for large space dimension d. 

The above conclusions are based on an analysis of the differential equations sat- 
isfied by the equal-time correlation functions of the scalar T. These equations may 
be obtained in many ways, for example with help of the MSR formalism. Let us first 
consider the example of the 2-point function. We have 

{T{t,Xi)T{t,X2)) = J (^J T{t,Xi)T{t,X2) e-^'-^'^'"') D R DT j noim^ d^i^{v) 

where 

S{R,T,v) = ij R{dtT + V VT - KAT)dtdx + \jR{t,x)C(^)R{t,y)dtdxdy 

and dn^lv) is the Gaussian measure corresponding to the covariance (122|). The normal- 
ization of the Gaussian R, T functional integral may be argued to be f -independent. 
Performing this integral, one obtains 

(T(t,xi)T(t,X2)) 

= J [{dt + V ■ V - KA)-^C{~dt + V ■ V - kA)-^) {t,xi;t,X2) d^i^{v). 
Expanding into the Neuman series 

oo 

{±dt + v-V - kA)-i = (±^t - ^^^y^ (-^ ■ V {±dt - kAY^Y , (28) 

m=0 

we may easily compute the rfyU^ {v) expectation using the Wick theorem. The resulting 
sum may be identified|^ with the Neuman series for M.2^C where M.2 is an operator 
acting on the functions of Xi,X2, 

M2 = -«:(A,,+A,J - ^(((^(t,a;i)- V,j2^ + ((i;(t,X2)- V,.J2)) 

- Ht, xi) ■ V^, v{t, X2) ■ v^.^) 

= -{^ + Do)(A,, + A,.J - 2{D^5'^ - d'^ix, - X2))d,.d^,^ . 

The Neuman series develops M.2^ powers of the w-expectations. The "tadpole" term 
{{v{t,Xi) ■ VxJ^) originates from the Wick contraction of 2 velocities in the Neuman 
series of the same {±dt + f ■ V — kA)~^ operator whereas {v{t, xi) ■ V^i v(t, X2) ■ V^a) 
comes from the mixed contractions. 

Higher point correlation functions can be treated the same way: 

J^2n{Xi, . . . ,X2n) = {T {t, Xi) ■ ■ ■ T (t , X2n) ) 
» n 

= E y n((^* + ^ ■ V - i^^y'ci-dt +V-V- ka)-!) (t, x^r, t, dfi^{v) 

pairings j — \ 
({"j >™j }) 



and using eq. (pH]) and the Wick theorem, it is easy to obtain the inductive relation 

J^2niXi,...,X2n) = J2 ^ . _ . . X2n) (29) 

l<ni<n2<2n "a 
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the time-decorrelation of the velocities is crucial at this point 



(in a slightly abusive notation) with 



J 

m=l l<m\<m2<n 



For K > 0, /An are positive 2"^^^ order elliptic operators. For k = they become 
singular elliptic (their principal symbol looses positive definiteness at coinciding points 
= Xm2)- Note that in the action on translationally invariant functions of Xi, . . . , a;„, 

n 

Mn = -nYl Ax™ + 2 E ^''(^-1 - . (30) 

m=l l<»ni<m2<n 

Hence in the translationally invariant sector the constant Dq ~ divergent as the 
IR regulator m —>■ decouples from operators For k, = and m = 0, turn 

into scaling operators Ai^^ of homogeneity degree ^ — 2. 

The 2-point function equation 

•F2(xi,X2) =A^-iC(^) (31) 

may be rewritten (with the scaling form of d^'' and r = |si — a;2| ) as 

(^2^2)(r) = 2ri-'^9, (Z^/-^+« + ft:/-^)9,^2(r) = C(9 (32) 
leading to the explicit solution 

~ JC{l)a''-^da 
^^W = iy + -iP. (33) 

r 

The integration constants were chosen so that drJ-'2{0) = and ^^2(00) = which is 
equivalent to the use of the Green function kernel for Ai2^ in eq. (0). Such a choice 
describes the stationary 2-point function obtained from a localized initial distribution 
of T by waiting long enough. It is trivial to analyze the integral in eq. ([33| ) explicitly 
and to see that, first, the k — ^ limit of JF2 exists and, second, that 

Mr) = A,L'-^-^r^-^ + O(L-) . (34) 

L large 

In particular, for the 2 structure function of T, we obtain 

S,{r) = 2^(0)-2^2(r) = i^i^r^'^ + 0{L-') 

in agreement with the naive dimensional analysis of eq. (0). By similar dimensional 
analysis of eq. (|35|) , we may expect that 5*2,1 (r) ~ r"^^~'^) for large L. This is the 
Kolmogorov-type prediction and it is incorrect for n > 1. The first hint of what may 
go wrong with the dimensional argument may be already seen in eq. (^4|) where on 
the right hand side, beside the scaling contribution — 2Dd{2-i) ' ^^ere appears the 
constant AqL'^"^ diverging when L oo. Any divergent contributions to J^2 have to 
be annihilated by A^2 since the right hand side of eq. (0) is regular when L oo and 
clearly constants are zero modes of Al2- 

Similarly, contributions annihilated by A^2n may appear in the (2n)-point correla- 
tion function JF2„. A more refined analysis shows that, for sufficiently small ^ and for 



fixed Xi, . . . , X2n, the k — > and m — > limits of exist and are dominated for large 
L by the contributions of the scaling zero modes of Aif^'- 

^2n(x) ^ = ^, Ac,.„^"^'-^^-^^"^2n(x) + + . . . . (35) 

L large 

where jFg^ is a zero mode of J^fn of the homogeneity degree (2n given by eq. (^7^ and 
the dots denote terms which do not depend on one of the variables Xm and consequently 
do not contribute to the structure function S2n- The amplitudes Ac,2n are non- universal 
in the sense that they depend on the shape of the forcing covariance C. The zero mode 

•^2n(x) = E (3^1 - 2:2)' ■ ■ ■ ix2n-l " X2n? + 0(0 + • • • • 

permutations 

and reduces for ^ = to a polynomial zero mode of the (2n(i) -dimensional Laplacian. 
The essential tool in arriving at the result (P5|) is the use of the Mellin transform of 



eq. (^) with the scaling operator (A^2n)^^, supplemented by the perturbative expan- 
sion of the scaling zero modes of A4^j^ in powers of ^. The analysis has a renormalization 
group flavor. The perturbative argument is applied to the single scale problem in which 
the differential operator restricted to scaling functions of a given homogeneity de- 
gree is analyzed. Such an operator has a discrete spectrum. The perturbative zero 
mode information is then plugged into the inverse Mellin transform which assembles 
different homogeneity degrees. More exactly, the scaling zero modes enter the residues 
of poles of the Mellin transform of the Green function (A^2n)~^- The above analysis 
should be contrasted with the direct perturbative expansion of J^2n in powers of ^ which 
requires perturbative treatment of the Green function of A^l*^, an operator with a con- 
tinuous spectrum in L^. The latter expansion is plagued by logarithmic divergences 
(proportional to powers of logL) which are resummed on the right hand side of (|35|). 
Still, the anomalous 0{^) contribution to (2n may be extracted from the coefficient of 
logL in the 0{^) term of J-'2n proportional to the integral 

„ iip-x+k-y+q-z)/L (p-k){q-k) 

J P' 1' p2 _^ J^2 _^ q2 _ . _^ ^^dp2q2 

where Sp^kfip) = f{p+k)—f{p). We shall unravel the renormalization group underlying 
the exponentiation of the logarithmic divergences in the perturbative treatment in 
powers of ^ of the passive scalar in the next section. 

The correlation functions J^2n are not smooth at coinciding points even for k, > 0. 
Nevertheless 

lim {VT){x) ■ {VT){y) = e{x) 
y 

exists inside the correlations for k > and defines the dissipation field. In particular, 
the mean dissipation rate e = (e(a;)) = 2^(0), as may be easily seen from eq. (|33|). 
t{x) does not vanish when /t ^ but is given by the dissipative anomaly 

lim e(x) = lim d^{x - y){^.^^T){x) {dyjT){y) (36) 
K — i-0 y — *x 

holding inside correlation functions. The result (B^) permits to infer the scaling behav- 



ior ( p6| , p7D of the structure functions and, together with the dissipative anomaly (|36|), 
also the inertial range scaling of the correlations involving the dissipation field. One 
obtains, for example. 



(e(x)e(y)) -e^ ~ ^4-2(2-0 



LECTURE 4 



It has been realized a long time ago |31|] that there exists a similarity between the 
behavior of the 2-point correlation functions of a nearly critical statistical-mechanical 
systems and of the Fourier transform of the equal-time velocity correlators in a turbulent 
flow. Both have a scaling regime of power-law decay followed by much stronger decay 
at infinity. One may then establish the following dictionary ||35|in 



critical phenomena turbulence 

UV cutoff integral scale 

inverse correlation length viscous scale 

T — Tc viscosity v 

scaling regime inertial range 

anomalous conservation laws dissipative anomaly 



This suggests that, very roughly, the turbulant phenomena look like critical phenomena, 
provided that we invert the scales interchanging short and long distances or the position 
and the Fourier spaces. Were this true, the short- distance universality of the critical 
phenomena (independence of the long-distance behavior of the microscopic details of 
the system) should be accompanied by the long-distance universality in turbulence 
(insensitivity of the short- distance behavior to boundary effects or/and details of the 
energy injection). 

The right tool to study the scaling properties of the critical phenomena and to 
establish their short-distance universality has been provided by the Kadanoff- Wilson 
renormalization group (RG) [^, Loosely speaking, the RG idea is to look at the 
system from further and further away so that its microscopic details are eventually 
wiped out and many microscopically different models start looking the same. By anal- 
ogy, it seems |3^, ^ that the turbulent systems require an inverse renormalization 
group (IRG) analysis. By examining them through a stronger and stronger magnifying 
glass, we would loose their large-scale details from the vision range and should discover 
a short- distance similarity of different turbulent cascades. The presence of a finite but 
large correlation length in the nearly critical systems beyond which there is a crossover 
to the high temperature regime would then correspond to the presence of the short 
viscous scale in high Reynolds number flows beyond which the dissipative regime sets 
in. 



RG had an enormous success in explaining critical phenomena |^T], Why is 
it then that IRG never developed beyond the level of a vague idea? Is the analogy 
between the critical phenomena and turbulence too naive and missing totally the es- 
sential points? In the author's opinion the reason is different. RG is not a universal key 
to all problems as it is sometimes thought. Its effective use requires a correct choice 
of RG transformations and that, in turn, requires a good understanding of physics 
of the system. Similarly, an IRG-type analysis of turbulence will require a deep use 
of knowledge of physics of turbulence. To provide an argument for such a thesis, we 
shall show that the IRG idea allows to systematize the analysis of Kraichnan's passive 
scalar described in the previous lecture and opens a possibility to extend it to more 
complicated systems. 

^^see also |l3| for a more recent discussion and more references 



Let us start by a short reminder of how one may perform a RG analysis of the 
long-time long-distance asymptotics of the (nearly) critical dynamics described by the 
stochastic evolution equation 

dtT = DAT - ^ + / (37) 

with V(T) = JV{T(t))dt and the Gaussian noise / with mean zero and covariance 

{fit,x)f{s,y)) = S{t-s)C{x-y). 
The corresponding MSR action is 

S{R,T) = So{R,T) + Si{R,T) 

where 

So{R,T) =ijR{t,x){dtT -uAT){t,x)dtdx + \j R{t,x)C{x - y)R{t,y)dtdxdy, 
S,{R,T)=rjR{t^x)^^dtdx + \ j J^^^^'^^ dtdx, 

see eq. (0). The Sq part of the action may be used to define the Gaussian "measure" 

dn^{R,T) = e-^o(^'^) Di^DTy/ norm. 

with covariance G given by the 2-point functions of the forced diffusion (|14|,|l5|). 

Let R, T and R, T be two copies of Gaussian fields distributed with measure dfic- 
We shall decompose 

R = Ri/x + P, 

T = Ti/A + r (38) 

demanding that p, r be Gaussian fields independent of R, T. The rescalings of R, T 
are as in eq. (16), i.e. they serve to exhibit the long-time long-distance scaling of the 
linear forced diffusion. By assumption, we have the factorization 

dp^{R,T) = dp^{R,T) dfi^^{p,T) . (39) 

The 2-point functions of p, r building the covariance Fa are the differences of the 2- 
point functions oi R,T and Ri/x, Ti/x, very much in the spirit of the Pauli-Villars 
regularization. Explicitly, 

{T{t„x,)T{t„x,)) = J el*-*^l«^^-Hxi-..) £^^Sk. (40) 

i.e. it is the high-momentum part of the covariance (TT). The 2-point functions 
involving p vanish (recall that {RR) = and (TR) is scale-invariant) but we shall 
keep p in the formulae which will be later applied in situations with p 7^ 0. The 
decomposition (RH) of T is into the low- momentum part Ti/x and the high- momentum 
fluctuation r and allows to define the effective interactions 



J Q-Sl{R,j,+P,T,j,+r) dp^^{p,T) (41) 



by integrating out the high- momentum fluctuations from the Boltzmann factor e 
The RG transformations 7l\ : Si S\ have a semigroup property, 7l\ o 71^' = T^w ■ 



We may also integrate out the high- momentum fluctuations in the insertions Fi [R, T) 
into the MSR functional integral, defining the effective insertions by 

F^{R,f) = I F^{Ry, + p,fy, + T) e-^i(^VA+P.?iA+-) ci;,^^(p,r)/e-^^(^'^). (42) 

Note that the transformation Fi i-^ Fx is really the derivative dTZx{Si) of the semigroup 
Tlx and that the RG transformations preserve the expectation values: 

{Fi)s, = (43) 
where (F)^^ = / Fe'^^ / J e"'^^ . 

The vanishing interaction Si = corresponds to a (Gaussian, trivial) fixed point of 
the semigroup Tlx- For non-trivial Si, the simplest situation occurs if under the action 
of Tlx it converges to a (possibly trivial) fixed point: 

Sx ^ — > . 

X—i-oo 

The study of the long-distance asymptotics of the correlation functions reduces then 
to the search of the corresponding scaling fields. Let us explain the latter concept. Let 
F{R, T; x) be a functional of R, T explicitly dependent on a sequence of space-points 
X. For example, we may take F(i?, T;x) = Y[T{t,Xi). Suppose further, that for some 

i 

exponent 

A-«*(F(Ax))a F,(x). (44) 

Then, tautologically, 

(F,(x))a = A«*F,(A-ix) 

if in the computation of the effective insertion on the left hand side we use the fixed 
point interaction 5*. Hence the name: scaling field for -F*(x). Besides, in view of eq. 

lim A-^*(F(Ax)),^ = (F.(x)),, 

A— >oo ^ 

giving the long-distance asymptotics of (F(x))g^ (if the fixed-point expectation on the 
right hand side does not vanish). 

The first information about the RG fiow Tlx Hiay be obtained by studying its lin- 
earization around the Gaussian fixed point. In fact for functionals 5*1 polynomially 
dependent on R, T with local scaling kernels 

Sx = dTlxiO)Si + 0{S^) = At^^l Si + lower order polyn. + 0{S^) (45) 

where the (long distance) dimension [Si] of Si is calculated additively with the use of 
the following table 

H = l, M = 2, [T] = l-I, [R] = -l-l. 

Small irrelevant interactions with [Si] < should then die out under the iterated RG 
transformations]^ resulting in the convergence of Sx to the trivial fixed point and in the 

^^this may require fine tuning of the lower order terms in Si 



long-time long-distance asymptotics of the correlations as in the linear forced diffusion 
case. The fate of relevant {[Si] > 0) or marginal ([5*1] = 0) interactions cannot be 
determined by the linear analysis around the trivial fixed point and requires higher 
order calculations which may show convergence to a non-trivial fixed point situated 
in a perturbative neighborhood of the trivial one. For example, the nonlinearity with 
V(T) ~ / T'^ in eq. (|37|) , describing the Langevin-type dynamics for the (j)'^ field theory, 



leads to the interaction ~ / RT^ of dimension — 1 — 1+3 — Y+2+(i = 4— (i, in agreement 
with the well known static power counting rendering the nonlinearity irrelevant above 
4 dimensions. For d < 4, the RG flow is instead governed by a non-trivial fixed point 



which seems accessible by a perturbative e-expansion in powers oi e = A — d |^2[. An 
important aspect, essential for the validity of the RG analysis in the above systems, 
is the approximate locality of the effective interactions Sx in the position space which 
means physically that no low-energy interacting modes were removed from the system 
by integrating out the short-distance fluctuations. Such locality (usually checked only 
perturbatively) allows to separate the dominant exactly local scaling contributions to 
driving the RG dynamics from the remainder strongly damped under the RG flow. 
Technically, the separation is done by Taylor-expanding the kernels in Sx in the Fourier 
space (the approximate locality of the kernels in the position space makes their Fourier 
transforms smooth). 

Suppose that instead of being interested in the long-time, long-distance behav- 
ior of the nonlinear forced diffusion, we want to study its short-time, short- distance 
asymptotics. In particular, we would like to know how the nonlinearity effects the UV 
asymptotics ([T9|) of the solutions of the linear equation. In order to study this problem, 
we may set up an (inverse) RG formalism in full analogy with the one described above, 
with the only difference that we shall use the field rescalings (|18D suitable for tracing 
the short-time, short- distance asymptotics of the linear forced diffusion instead of (p!6|) 
appropriate for the long-time long-distance behavior. Repeating the decomposition 



(38) leading to the factorization (|39|) for the new rescaling (marked in the notation by 



superscripts rather then subscripts), we shall obtain the fluctuation r covariance 



Note that C{k) — X~'^C{k/X) is the Fourier transform of the C{x) — C{Xx) so that it 
corresponds to the long distance part of C. In other words, the decomposition T = 
T^/^+r is now into the short-distance part and the long-distance fluctuation. Repeating 
the deflnitions (^) and (|i2|) for the new scaling, we obtain the IRG semigroup TZ'^ : 
Si = S^ ^ S^ with the derivative dn^{S^) : Fi = F^. If the effective long- 

distance interactions S'^ converge to a flxed point S*, the study of the short- distance 
asymptotics of the correlation functions for the perturbed version of the forced diffusion 
reduces to the search for the scaling flelds 

lim A^* (F(x/A))^ = F*(x) 

A— ^00 

for which 

hm A^*(F(x/A))^, = (F*(x)),.. 

For small interactions S^, 

= dn^{0)S^ + 0{{S^f) = A^t^'l + lower order polyn. + OUS^) (47) 



where the (short distance) dimension [S^] is calculated with the use of the new table 



[a;] = l, [t]=2, [T] = l, [R] = -l-d. 

Note the change of the sign in the exponent of A in eq. p7| ) as compared to eq. (^)). 
For example, nonlinearity V(T) ~ /(T)^ in eq. (|37|) leading to the interaction ~ / RT^ 
of dimension —1 — d + 3 + 2 + d = 4, i.e. irrelevant by power counting. / RT^ couples, 
however, to the unstable constant mode of T. Considering instead the gradient-type 
nonlinearity with V(T) ~ /(VT)^ in eq. (0), avoiding coupling to the constant mode, 
one obtains the interaction ~ / i?V(VT)^ of dimension —1 — d— 1 + 2 + d = 0, 
i.e. marginal in all dimensions. Hence the linearized IRG does not provide any simple 
hints about the short-time, short-distance asymptotics of the nonlinear forced diffusion. 
Besides, one should check that the IRG effective interactions 5''^ possess in this case a 
Fourier space locality properties which would allow to separate a finite number of scaling 
contributions driving the IRG flow. It should be also noticed that the covariance (|46|) , 
unlike its RG counterpart (^), is not positive which may lead to non-perturbative 
complications in stabilizing the IRG flow. We shall have nothing more to say about 
the UV regime of the forced linear diffusion except repeating that its control is an 
interesting open problem with physical relevance. Below, we shall apply the IRG idea 
to the passive scalar model of Lecture 3 with milder nonlinearities and milder stability 
problems. 

The mutual relations of RG and IRG lead often to a confusion stemming from the 
fact that the RG is also used to study the short-distance asymptotics in field theories 
governed by UV fixed points. Thus there are two contexts in which we apply the 
standard RG: either we fix the UV cutoff and study the long-distance behavior of the 
theory by observing stabilization of the system under RG which lowers the momentum 
cutoff (the statistical-mechanical context) or we start with theories with a larger and 
larger momentum cutoff and apply RG to lower it to a fixed value, adjusting the 
parameters of the cutoff theories as to obtain stabilization of the effective theories on 
the fixed scale (the continuum limit or field theory context). As is well known, the 
two contexts differ essentially only by a straightforward rescaling and RG used in both 
of them integrates out the degrees of freedom corresponding to the shortest distances 
present. Similarly, IRG may be applied to systems forced on long distances in two 
contexts differing essentially by a simple rescaling. Either we fix the size of the system 
and the forcing scale (the IR cutoffs) and study the short- distance behavior by trying 
to exhibit stabilization under IRG which lowers the long-distance cutoffs (the more 
common situation) or we consider larger and larger systems forced at longer and longer 
distances and apply IRG to lower the distance cutoff to a fixed value hoping to see 
stabilization of the effective theories on the fixed distance scale (the infinite-volume 
limit context). Again, in both contexts the IRG transformations lower the distance 
cuttoff. In a given situation whether RG or IRG should be used depends on which 
leads to a finite number of expanding directions in the interaction space^. It seems 
that the limitation of the energy injection to long distances condemns the RG analyses 
of turbulence employing a modified forcing with a power spectrum and the e-expansion 
techniques O, O, p[. 



We would like to study via IRG the short-distance asymptotics of the passive scalar 
l]) with the random velocity distributed according to eqs. ( ^2] - ^ ). Under the rescaling 



v\t,x) = X^-'v{t/X^-^,x/X), (4J 
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it is a priori not excluded that there are exotic systems which yield to both RG and IRG analysis 



for large A. Hence (vv) exhibits a scaling behavior at short distances, modulo a diver- 
gent constant zero mode. One of the problems with the passive scalar is that the term 
— V- VT couples T to the unstable constant mode of v. A possible solution is to replace 
eq. (^) by its quasi-Lagrangian version 

dtT = -{v -vo) -VT + kAT + f . (49) 

where vo{t, x) = v{t, 0). This changes the stochastic dynamics. The new djTiamics may 
be interpreted as describing the system in the frame moving with one of the material 
points of the flow ||^ . Although eq. (^) corresponds to different dynamical correlations 
of T, the equal-time correlations do not change |]18|. As in Lecture 3, one may obtain 
for the latter the same equations (^) with directly in the translation-invariant 
form (^). Of course, instead of employing the new equation (^), we may keep the 
original one (pT]) changing only the velocity covariance to 



{v\t,x)v\s,y)) = -26{t- s)[d'\x-y)-d'^{x)-d'^{y)) . (50) 
It will be more convenient to choose another (vv) covariance, which we shall call mixed, 

{v\t,x)v^{s,y)) = -6{t~s)(2d'^{x-y)-d'^{x)-d'^{y)) (51) 



which leads to the equations ( pof ) with A4n replaced by 

n 

m=l l<mi<m2<n 

Ai'^ again coincide with Ain in the translational-invariant sector so that the mixed 
distribution ( pi]) of v leads to the same equal-time correlators of T. The covariance 
( 1^ is 0{C,) for K = and its use will simplify the perturbative expansion in powers of 
^. Although it lacks positivity, unlike the original covariance or the quasi-Lagrangian 
one, this will not cause stability problems. 



We shall introduce another modification of eq. (|2T|) by replacing it by 



dtT = -vVT + kAT - r]m^T + f (52) 



where m is the same IR regulator that in eq. ( pi)) and r/ > is a small fixed constant. 
The addition of the "mass term" ~ m^T, which changes the long-distance behavior is 
innocuous for the short- distance one, see below. The expectation values for the passive 
scalar (|52|) with the velocity distribution (|5TD may be generated in the MSR formalism 
in the following way. Introduce operator 

Mm = -d'^{x)d^.d^j + r]m^ . 

Mm is a generator of inhomogeneous superdiffusion (on distance scales <^ m~^) and 
^-tMm describes the dynamics of the 1-point function of the passive scalar with velocity 
covariance (O). Let 



S°{R,T) = i J R{t,x){dt- Mm)T{t,x)dtdx + \ j R{t,x) C{x - y) R{t,y) dtdxdy, 



S^{R,T,v) = *y R{t, x)v(t, x) ■ VT{t, x) dtdx , 
S'^{R,T,v) = —in J R{t, x) AT{t, x) dtdx , 
S"\R,T) = I J R{t,x)MmT{t,x)dtdx. 



Let d/j, I (f ) denote the Gaussian measure with the mixed covariance (0). The equal 
time (2?7.)-point function jF2„(x) of the scalar T may be represented as the MSR func- 
tional integral 

•^2„(x) = / \\T{t,Xi) e-^^^'^''') DRDTdfi. (v) / norm. (53) 
i=i / 

where S = + + S''^ + S"^. We have included the —i J RM^T term into free 
action compensating it by the S"^ term treated as an interaction. The role of S"'^ 
is to remove the "tadpole" contractions {{v ■ V)^) in {{S^)^). The latter are instead 
resummed into the free (TR) = i{dt — Mm)'^^ propagators originating from 5*°. Reading 
the field short- distance dimensions from and, for v, from the short-distance scaling 
of the (vv) covariance, we obtain 

[x] = l, [t] = 2-e, [T] = l-ie, [R] = -l-d+'-^, H=e-1. (54) 

Note that the Kolmogorov value of the velocity exponent ^ = | is obtained by equating 
the dimensions of T and v which is encouraging in view of the fact that the v ■ VT term 
of the passive scalar equation is replaced by v ■ Vf in the NS equation. The table (^) 
gives 

[S'] = 0, [S"']=0, [s"] = ~c 



We may attempt an IRG analysis of the short-distance behavior of the correla- 
tions J-'2n, basing the IRG transformations on the Gaussian measure dfi^{R,T,v) ~ 
e~^°^^''^^ DRDT dfi^,{v) and rescalings corresponding to the table of dimensions (p^, 

R\t,x) = \-^-''^^'^R{t/\^-^,x/\), 

T\t,x) = \^-^'^T{t/\^-^,x/\) (55) 

and v^it^x) as in eq. (|48|) . We could hope that the effective interactions obtained 
from the marginal terms + S"'^ tend to a fixed point, with the relevant S' term with 
very small n destabilizing the IRG trajectory only at very short distances causing 
eventually a crossover to the dissipative regime. A closer analysis of the effective 
interactions based on the ^ expansion shows however a lack of convergence of defined 
this way to a fixed point [Q. It appears that an infinite number of relevant terms is 
generated which destabilize the effective trajectory. In such situation, it may seem that 
IRG fails to predict the short-distance scaling of the scalar correlations in the inertial 
range. How is it possible then that we still were able to control this scaling for small 
^, as discussed in Lecture 3? 

There appears to exist a simple solution to the above paradox. The idea is to exclude 
the ^ / RCR term from the free action expanding in the functional integral (Q) 

e-\!^R= f with v4„ = ^ and FJR) = JR{t,x)C{x-y)R{t,y)dtdxdy. 

n=0 

Note that only the n^^ term of the expansion gives a non-zero contribution to (^) 



(the number of the T and R insertions must be equal now). The leftover free action 
S"° = i J R{dt — Mm)T leads to the Gaussian measure with the 2-point functions 
{RR) = = {TT) and {T{ti)R{t2)) = i{dt- Mm)-^{ti,t2) = - ta) e-^^^-^^)*^™ . It 
will be more convenient to introduce a cutoff version ci/x ^ (i?, T) of this measure with 
the 2-point functions 

{RR) = 0, (TT) = 0, 

(r(ti,xi) R(t2,X2)) = i{dt + Mm)~^(ti,xi; ^2,^2) - i{dt + M^rnY^iti.xi] 12^X2) 

= iOih - t2) (e-^*!-*^)^^'" - e-(*i-*2)^'^-) (xi,X2) . 

For A ^ 00, {T{ti)R{t2)) tends toQ - 12) e"(*i"*2)M„ except for ti = ^2 for which 
it always vanishes. Finite A introduces a short-time, short-distance cutoff into the 
heat kernel e"*^*^"*^-'^^'" . The functional- integral expression (|53|) for JF2„ may be now 
rewritten as 

2n 

/ n Tit,x,) F^iRT e-^s^+s'^)^R^T,.) d^^^,^{R,T) dfi^,Jv) 

J^2n(x) = An lim 

A^oo norm. 

2n 

^ A. (nT(t,x.)F,(i?r)^^^^,,. (56) 

i=l 

The role of the short-distance cutoff A is the same as that of the 5*"^ term before. It 
excludes the tadpole contractions {{v(t) ■ V)^) in {{S^)^) now forbidden because the 
cutoff {T{ti)R(t2)) propagator, unlike its A = 00 version, vanishes at equal times for 
all A. Note how eq. ( |5BD works for the 2-point function. Expanding on the right hand 
side e"'^^^'^' into the power series, computing the Gaussian expectations and passing 
to the A — > 00 limit, one obtains the Neuman series for Ai2~'^C resulting from treating 
as a perturbation the second line in the identity 

M'2 = {Mm)^, + (M„),, 

+ {2d''{xi - X2) - d'^ixi) - d''ix2))d,.d^,^ - kA,, - kA,, . 

The higher-point function formula works similarly. 

In order to generate the IRG transformations, we shall follow the general rules 
described before. Introducing the Gaussian fields R, T, v distributed with the measure 
dfi {R, T) djj, , iv) and decomposing 

R = R^/^ + p, T = f^'^ + T, v = v^'^ + w 
and, for the measures, 

d^^^A (R^ T) = dfi {R, f) dfi (p, r) , dfi (v) = dfi (i!) dfi (w) 

with 

= n' — n' 

^^the rim? term in M„ plays here the crucial role 



we shall define effective interactions 5''^ = TZ-^IS^) by integrating the long distance 
fluctuations p, r, w in the Gibbs factor e"*^ , 

The effective insertions will be given by the linearization dTZ'^ of the above IRG flow 
for interactions, compare eqs. ( p|^2D . Denote {T{t,xi) — T(t,X2))^" = F2„(T;x). 

Our main claims are as follows: 

1. For the molecular diffusivity k set to zero, the effective interactions tend to the 
(non-trivial) limit 5* as A — > oo, at least order by order in the expansion in powers of 
S^. At a fixed order, S*{R, T, v) is given by an explicit integral expressions with kernels 
approximately local (i.e. fast decaying) in the Fourier space. 

2. The limit 

lim A"(2-« (F2„(x/A))^ = F;„(x) 

exists and is a scaling field of dimension n{2 — ^) of the fixed-point theory. 

3. The limit 

lim A^2"-"(2-«) (F")^ = (F?)* 
exists and is a scaling field of dimension (2n — n{2 — C) = ~ '^^d+2^ ^ ^(C^); see eq. 

4. Finally, the limit 

^lim A^- (F2„(x/A)FY = F:^(x) 
exists and is a scaling field of dimension C,2n- 

The last three results have been established to the first two non-trivial orders in ^. We 
shall not discuss here the details of the perturbative analysis referring an interested 
reader to 0. Let us only mention the role of the Mellin transform of the kernels 
entering the effective interactions or insertions in separating the contributions with the 
lowest dimensions from the reminders strongly damped by the IRG flow. 

What are the implication of the above results? First, note that, in view of eq. (|56|) , 
the relation 

hm AL(F2„(x/A)F,)^, = (F^:,(x)),. (57) 

implies the anomalous scaling (|26|) : unlike the expectations {F2n) s* ((-^")*)s* 
which involve unequal numbers of R and T flelds, the right hand side of eq. (^7|) does 
not vanish. The even more interesting observation is that the anomalous dimensions 
C2n — n{2 — ^) are carried by the scaling flelds {F^)* which are relevant (i.e. of negative 
dimension) for small ^ and which correspond to the composite operators F". The 
scaling zero modes of the operators M.2n which played the crucial role in obtaining (p6|) 
enter the kernels in the explicit formulae for (F")*. Recall that in statistical mechanics 
or fleld theory local composite operators are produced by multiplying flelds localized 
at the same space-point. In the spirit of the scale inversion discussed at the beginning 
of this lecture, composite flelds in turbulence should be obtained by multiplying flelds 



localized at the same wavenumber in the Fourier space. Since the forcing covariance 
C is almost a constant in the inertial ranee, ^ C(0) / R{t, k)\^^dt so that the 
insertions are almost local in the wavenumber space. On the other hand, the scaling 
fields -F2*„(x) corresponding to the insertions (T(i:, Xi)— T(t, X2))^"' carry the Kolmogorov 
(normal) part of the dimension of the structure functions. There are no extra anomalous 
dimensions appearing in the scaling fields F * ^ (x) which correspond to the products of 
F2„(x) and F". Their dimensions are just the sum of the dimensions of F2*„(x) and of 
the composite scaling fields The IRG picture of the system permits to analyze 

systematically such operator products in the spirit of operator product expansions of 
the long distance type (as contrasted with the field-theoretic short-distance OPE's 



resulting from the RG analysis). Such long-distance OPE's lead to fusion rules of the 



type studied in [28 



Up to now, we have ignored the relevant 5"^ contribution to the MSR action of the 
model, proportional to the molecular diffusivity k. As discussed above, it eventually 
causes a crossover of the IRG trajectory from the vicinity of the convective fixed point S* 
to another one corresponding to the dissipative regime dominating very short distances. 
It should be interesting to study the crossover and the dissipative regime using IRG. 
As for the role of the rjm'^T term which we have added in the passive scalar equation 
(1^, note that its change gives rise to a term ~ i / RT in the MSR action of dimension 
2 — i.e. irrelevant at short distances. We could, in fact, analyze directly the 1] = 
case, introducing rj > only in the intermediate IRG steps. 

Summarizing. The IRG analysis explains the breakdown of the Kolmogorov scaling 
of the higher structure functions of Kraichnan's passive scalar for small ^ as due to the 
appearance of relevant composite scaling fields (-F?)* with a multifractal spectrum of 
dimensions. These fields do not destabilize the convergence of the effective interactions 

to the fixed point S* since their appearance in S'^ is forbidden by the conservation 
law imposing the equality of the numbers of R and T fields in interactions (this was 
not the case in our initial attempt to perform the IRG analysis which kept the J / RCR 
term in the action). Note the similarities but also the differences with the RG picture 
of the critical 0^ theory below 4 dimensions governed by a non-trivial 0{e = 4 — d) 
fixed point with a modified scaling^ and a finite number of relevant scaling fields. We 
believe that although our analysis has used crucially the simplifications inherent in the 
Kraichnan's model, the above conclusions are robust enough to explain the occurrence 
of multifractal exponents in other turbulent systems, for example in the passive scalar 
model with velocities correlated in time [P]. Whether a similar mechanism is responsible 
for the breakdown of the Kolmogorov scaling in the full-fiedged NS turbulence remains 
to be seen. One of the main problems there is the lack of a small parameter (like 
^ above) which will render obtaining reliable numerical values of structure-function 
exponents very difficult, if not impossible. The situation may be different in weak 
turbulence P3| which can provide a more gratifying test ground for the IRG ideas. 
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